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Mh The X(3872) is studied as an axial- vector charnionium state in the mul- 

O , tichannel framework of the Resonance-Spectrum-Expansion quark-meson 

D model, previously applied to a variety of other puzzling mesonic resonances. 

,S^ Included are the open-charm pseudoscalar-vector and vector-vector chan- 
nels, the most important of which is the S-wave D*^D^-\-D*^D^ channel, 

,__! which practically coincides with the X(3872) structure. The two free pa- 

^ rameters of the model are tuned so as to roughly reproduce the Xci(3511) 

\^ mass as well as the enhancement just above the I)*^D^ /D*^D^ threshold. 

OO The present model is able to describe the shape of the latter data quite 

■^ well. However, as no dynamical resonance pole is found, the X(3872) and 

^^ X(3940) cannot be reproduced simultaneously, at this stage. A possible 

l/~j further improvement is discussed. 

o 

O PACS numbers: 14.40.Pq, 13.25.Gv, 11.80. Gw, 11.55.Ds 



1. Introduction 

The X(3872) charmonium-like state was discovered in 2003 by Belle [1], 
as a -K^TT^ J/ip enhancement in the process B — )• K tt^ ■k~ J / if) . The same 
enhancement was then seen in several other experiments as well. Finally, 
the X(3872) was also observed in the Z)°L»°7r° and D*^D^ [3] channels. 
It is now listed in the PDG tables [4J, with a mass of 3872.3 it 0.8 MeV and a 
width of 3.0lti'4 ±0.9 MeV. The observed decay modes and helicity analysis 
indicate positive C-parity and favour 1++ or 2 ^ quantum numbers ^4]. 
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Theoretical work on the X(3872) comprises a variety of model calcula- 
tions and interpretations, because of its seeming incompatibility with the 
standard charmonium spectrum. Let us briefly summarise what we consider 
the most significant approaches. 

Tornqvist [5] Calls it a deuson, a deuteron-like DD* bound state. Binding 
gets provided by pion exchange. Large isospin breaking is expected, 
and so an important J/ip p decay mode 

Bugg ^ Examines a mechanism for resonance capture at thresholds that 
involves a threshold cusp interfering constructively with a resonance 
from confinement and/or attractive t and u-channel exchanges. Con- 
clusion: the X(3872) requires a virtual state or a resonance 

Gamermann & Oset [7J Dynamical generation with isospin breaking, in 
an effective chiral approach. Preference for a slightly unbound, virtual 
D*^D^+D*^D^ state. 

Lee et al. [8j Possible DD* molecular bound state in a potential model 
from heavy-hadron chiral perturbation theory. 

Kalashnikova & Nefediev [9j Analysis from data concludes preference 
for a dynamically generated virtual DD* state, with a sizable 2 ^Pi cc 
component. 

Danilkin & Simonov |10| Employs the Weinberg eigenvalue method for 
a model with 1 confinement and 1 meson-meson channel, coupled via 
a ^Pq pair-creation vertex. Two ^Pi peaks are found, i.e., a narrow 
one at 3.872 GeV and a much broader one at 3.935 GeV, a candidate 
for the X(3940). The X(3872) is of a dynamical nature here. 

Bugg [llj Discusses several approaches, including the covalent bond |llj . 
Proposes resonances as linear combinations of qq and meson-meson. 

Besides these works, there exists a plethora of exotic models, which never- 
theless are incapable of describing the X(3872) production data. 

2. Resonance-Spectrum Expansion 

The Resonance-Spectrum-Expansion (RSE) [12J coupled-channel model 
is designed to describe scattering processes of the form AB — )• CD. Al- 
though A, B, C, and D can in principle be any structures, we usually con- 
sider only non-exotic mesons. Between the initial and final meson-meson 
phase, there is an intermediate one that we call the confinement phase, and 
it is built up from an infinite series of qq states carrying the same — again 
non-exotic — quantum numbers. 
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To describe this picture, we define the effective potential 



Vij{pi,Pj;E) = X'^Jl^iPia)ji(pja)Y^ 



9i{n)gj 



n] 



n=0 



E — En 



(1) 



where pi and p'- are the relativistically defined relative momenta of initial 
channel i and final channel j, and j^. is the spherical Bessel function. As free 
parameters we have A, which is an overall coupling constant, and o, a param- 
eter mimicking the average string- breaking distance. The relative couplings 
gi{n) between meson- meson channels i and qq recurrences n, with discrete 
energies En, are computed with the formalism developed in Ref. [14J. The 
transitions mechanism is the creation/annihilation of qq pairs, according to 
the OZI rule, using a spherical Dirac-delta potential, which Fourier trans- 
forms into a spherical Bessel function in momentum space. Due to the 
separable form of potential (IT]), the corresponding T-matrix is obtained in 
closed form, via the Lippmann-Schwinger equation with relativistic kine- 
matics [T2]. Note that the resulting S'-matrix is manifestly unitary and 
analytic. 

In Eq. M, we take a harmonic-oscillator (HO) potential with constant 
frequency, not only because of its simplicity but also its success in many phe- 
nomenological applications (see Ref. ^5] for references). The corresponding 
energy levels are 

En = mg + nig + u{2n + 3/2 + I) . (2) 



3. X(3872) as an axial- vector charmonium state 

In the present study we assume the X(3872) to be the 2^Pi cc state, 
i.e., the first radial excitation of the Xci(3511) [3j, with axial-vector (l"*""*") 
quantum numbers. Then, the most important OZI-allowed decay channels 
are pseuscalar-vector and vector-vector, shown in Table [TJ Note that, for 
notational simplicity, we have omitted the bars over the anti-D mesons. 



channel 


relative L 


threshold [MeV] 


DsDt 







3872 
3880 
4080 


DsDl 


2 
2 
2 


3872 
3880 
4080 


D*D* 


2 


4018 



Table 1. Included meson-meson channels 
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4. Results and conclusions 

Before studying resonance poles in the vicinity of the X(3872), we must 
first fix the model parameters. For the charm quark mass rric and the 
universal HO frequency u, we take the values 1.562 GeV and 0.19 GeV, 
respectively, which have been kept unaltered in all previous work, starting in 
Ref. [T^. As for a, one of the two free parameters, we take 2.0 GeV~^ (about 
0.4 fm), which is roughly in agreement with the value chosen in Ref. [15] 
for ss vector resonances, scaled with an expected factor of y/rrujrric. The 
overall coupling A we use here as a tunable parameter in order to study 
the behaviour of complex-energy poles, which can be found numerically 
with ease, since the T-matrix is known analytically. An indication for the 
approximate value of A may come from the mass of the Xcii^P) m, viz. 
3511 MeV. 

Table [2] shows some pole positions for different values of A. The poles 



A [GeV-3/2] 


Pole 1 [MeV] 


Pole 2 [MeV] 


3.2 


3551 


3871 


3.0 


3555 


3873 - il 


2.2 


3572 


3884 - i6 


1.2 


3590 


3928 - il6 



Table 2. Pole positions as a function of A 

labelled 1 and 2 originate from the HO confinement states at Eq = 3599 MeV 
and El = 3979 MeV (cf. Eq. ^ with i = 1), respectively, for A ^ 0. 
These we call confinement poles. We also searched for dynamical — or 
continuum — poles, which cannot be simply linked to a confinement state 
and, moreover, disappear in the continuum for A— t-O, with an ever increasing 
negative imaginary part. Such poles we observed e.g. in Ref. [15], though 
with a large width. Here, if one showed up close to the lowest threshold, 
just as in Ref. [lOj, it might explain the X(3872) structure as a dynamical 
resonance, but none was found. We shall come back to this point below. 

Focusing now on Pole 2, which comes out close to the X(3872) struc- 
ture for A ~ 3.0 GeV~^'^, we plot in Fig. 1 its trajectory around the 





Re^ 


Im k 


Re k 


type of pole 


1 


>Th 


<0 


>0 


resonance 


2 


<Th 


<0 


>0 


virtual resonance 


3 


<Th 


<0 


= 


virtual bound state 


4 


<Th 


>0 


= 


bound state 



Table 3. Pole 2 in different Riemann sheets (see text and Fig. 1) 
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D*'^D°/D*^D° threshold, through 
two different Riemann sheets (RSs). 
The trajectory can be spht into 
four branches according to the clas- 
sification in Table [3[ where Th 
stands for the threshold energy. We 
change from one RS to another by 
choosing a specific root of the rela- 
tive momentum k. The first three 
branches lie on one RS, and the 
fourth on another RS. We denote 
branch 2 by "virtual resonance" , as 
it concerns a resonance with nega- 
tive phase space. This is a phenome- 
non unique to 5-wave thresholds. 
For other waves, a resonance pole 
always turns into a virtual bound 
state directly, without the virtual-resonance branch. 

In Fig. [2| left-hand graph, we depict the elastic D^D*^ amplitude just 
above threshold, for 4 values of A corresponding to the 4 situations described 
in Table p^ and in Fig. 1. As a matter of fact, we plot the quantity A;|Tp, 
which is the one to be compared [7] with experiment. The scarce data, 
not shown here, cannot really distinguish between the 4 scenarios described 
above, though we conclude that a bound state is less likely [16]. 



Fig. 1. Trajectory of Pole 2: from a 
resonance to a bound state. 
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Fig. 2. fc|T|2 for elastic D'^D*° scattering (0 MeV corresponds to 3872 MeV). 
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Finally, in Fig. [2j right-hand graph, we generate a cusp structure around 
the D D*^ threshold, by reducing A and so moving Pole 2 upwards. While 
quite illustrative, such a scenario is now excluded by the data. 

In conclusion, we have studied the X(3872) enhancement as the 2^Pi 
cc state in the RSE coupled-channel model. The included OZI-allowed 
channels are capable of reproducing the observed structure in D^D*^, via 
a nearby (virtual) resonance or virtual-bound-state pole. However, the 
present model does not generate any dynamical pole close to the real axis, so 
that the X(3940) [H — if indeed a 1++ state as the X(3872) — remains un- 
explained. We are now working [16J on model extensions, e.g. by including 
the probably observed [17], OZI- violating J/tp p decay mode. 
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